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Direct photons from Au+Au collisions at RHIC: QGP vs. hot hadronic gas 
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We have analysed the preliminary PHENIX data on the transverse momentum distribution of 
direct photons in 0-20% centrality Au-|-Au collisions at ^sjviv=200 GeV. In ideal hydrodynamics, 
data are explained if Au-|-Au collision produces Quark-Gluon-Plasma at the temperature Ti=AQQ 
MeV, at an initial time Ti=0.6 fm. PHENIX data are not explained in the alternate scenario when 
Au-|-Au collisions produces hot hadronic gas at temperature Ti < 220 MeV at an initial time n < 

5 fm. 
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In Au-|-Au collisions at RHIC, one observe a dramatic 
suppression of high px hadrons Hiili . The sup¬ 
pression is more in central than in peripheral collisions. 
It is also established that high px suppression is a final 
state effect, a parton before fragmenting in to hadrons 
suffers energy loss in a dense matter, leading to sup¬ 
pressed production [^. High px suppression together 
with the observation that bulk of the hadron produc¬ 
tion data {px spectrum and elliptic flow et(^ are well 
described in a ideal hydrodynamic model |6|, l7| , strongly 
support the idea that RHIC has produced thermalised 
matter at very high energy density. However it is not 
certain whether the matter produced is strongly inter¬ 
acting Quark-Gluon-Plasma (sQGP) as predicted in lat¬ 
tice QCD calculations. Hadrons, being strongly inter¬ 
acting, are emitted from the surface of the thermalised 
matter and carry information about the freeze-out sur¬ 
face only. They are unaware of the condition of the in¬ 
terior of the matter and can provide information about 
the deep interior only in an indirect way. In a hydro- 
dynamic model, one fixes the initial conditions of the 
fluid such that the ’’experimental” freeze-out surface is 
correctly reproduced. In contrast to hadrons, photons 
being weakly interacting, are emitted from whole volume 
of the matter. Throughout the evolution of the matter, 
photons are emitted. Conditions of the produced matter, 
at its deep interior, are better probed by the photons. 

Recently PHENIX collaboration published their mea¬ 
surements of direct photons in different centrality ranges 
of Au-I-Au collisions at yisNN=200 GeV 0. pQCD 
model calculation of direct photons in p-bp collisions, 
multiplied with the thickness function, correctly repro¬ 
duces the data in all the centrality ranges of collisions. 
Transverse momentum distribution of direct photons in 
different centrality ranges of collisions scale with the 
number of binary collisions. Apparently the direct pho¬ 
ton results are in direct contradiction with high px sup¬ 
pression observed in Au-|-Au collisions. Partons frag¬ 
menting in to photons do not suffer energy loss in con¬ 
trast to partons fragmenting in to hadrons. Moreover, 
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very good description of the published data with pQCD 
photons donot leave any room for thermal photons, which 
are expected to be emitted in large numbers from the 
thermalised matter. 

Measurement of direct photons is very challenging, 
more so at low px- The huge background from 7r° de¬ 
cay, dominate the spectra and proper algorithm for back¬ 
ground subtraction is very important. Very recently 
PHENIX collaboration improved upon the analysis tech¬ 
nique for photon measurement. Compared to the conven¬ 
tional method the new method improves both the sig¬ 
nal to background ratio and the energy resolution. In 
QM2005, PHENIX collaboration presented the (prelim¬ 
inary) result of their new analysis With the new 
analysis technique, direct photon yield in the low px 
range (1-4 GeV), in 0-20% centrality Au-|-Au collisions 
is increased substantially. pQGD predictions no longer 
can explain the data. Direct photons in excess of pQGD 
(hard) photons are possibly from a thermal source. The 
PHENIX (preliminary) direct photon data thus provides 
the first opportunity to measure the initial condition of 
the matter produced in RHIG Au-bAu collisions, at deep 
interior. 

In the present paper, in a hydrodynamic model, we 
have analysed the PHENIX (preliminary) direct photon 
data. Procedure for obtaining photon spectra in a hydro- 
dynamic evolution is well known ■ We have solved the 
hydrodynamic equations = 0 for a baryon free gas 

assuming cylindrical symmetry and boost-invariance. We 
have considered two possible scenarios: (i) Au-I-Au colli¬ 
sions produce QGP as the initial state and (ii) Au-bAu 
collisions produces hot hadronic gas. In the first sce¬ 
nario produced QGP expands, cools, undergoes 1st or¬ 
der phase transition at critical temperature (Tc), enters 
a mixed phase, remain in the mixed phase till all the 
quark matter is converted into a hadronic matter then 
cools to freeze-out temperature. In the second scenario 
the hot hadronic gas expands, cools till the freeze-out. In 
both the scenarios, photons are emitted throughout the 
evolution and their yield in integrated over the space- 
time volume. The second scenario is very important for 
unambiguous detection of QGP. Direct photons are very 
strange probe. Theoretical predictions indicate that 
they are emitted equally well from the QGP and from 
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the hot hadronic phase. Only in a narrow transverse mo¬ 
mentum window around 3 GeV, the two phases may be 
distinguished. Indeed, direct photons measured at SPS 
energy in S-|-Au and in Pb-|-Pb collisions mm also 
raised high hope of detecting QGP. However, it was later 
found that the SPS ener gy data are well explained in 
models without QGP [ 13 . 11 ^. 

In the present calculation we have used the bag model 
equation of state for the QGP phase, Pq = aqT'^ — B with 
Qq = 42.257r^/90. The hadronic equation of state was 
generalized to include all the mesonic resonances with 
mass < 2 GeV. The cut off 2 GeV is rather arbitrary and 
we verify that the results do not depend on the value of 
cut-off significantly. The bag constant B was obtained 
from the Gibbs condition pqgp{Tc) = Phad{Tc)- The 
critical temperature (Tc) and the freeze-out temperature 
(Tf) are assumed to be 180 MeV and 100 MeV respec¬ 
tively. 

Photon emission rate from QGP and from hot hadronic 
gas are well known. For the single photons from hadronic 
gas we include the following processes, 

(a) TTTT —>■ p7, (b) np —>■ Try, (c) u —>■ Try, (d) p TTTry 
(e) np ^ Ai ^ Try 

rates for which are calculated in [lallJ. Photon emission 
rates from QGP are calculated in |l SL ll 3 mm 

Hydrodynamic models require initial time (rQ and ini¬ 
tial energy density profile {si{r)). n is essentially the 
thermalisation time beyond which hydrodynamic is ap¬ 
plicable. RHIC data on elliptic flow as well as hadron 
Pt spectra require ve^ short time scale of thermali¬ 
sation Ti=0.6 fm [ 3 . In the first scenario, when 
QGP is produced in the initial state, we use this value 
as the initial time. The PHENIX (preliminary) direct 
photon data are in the 0-20% centrality range of col¬ 
lisions. In this centrality range, the average number 
of binary collisions is < Nunary >~ 779. We have 
assumed that 0-20% centrality Au-|-Au collision corre¬ 
sponds to central collisions of nuclei Agff, such that the 
average number of binary collisions is reproduced. With 
TAAib = 0)- :^ = we obtain A^ff w 143. 

For the initial energy density profile we use a Woods- 
Saxon form with radius i?e//=5.85 fm (corresponding to 
nucleus Ae//=143). For the diffuseness parameter we 
use a = 0.54 fm. The central value of the energy den¬ 
sity could be obtained by fitting the PHENIX prelimi¬ 
nary data. However, presently we donot attempt exact 
fit to the data. Rather we have calculated the thermal 
photon yield with three values of central energy density, 
e(r = 0)=27.8,47.1 and 75.1 GeV/fm^ corresponding to 
initial temperature Ti=0.35, 0.40 and 0.45 GeV respec¬ 
tively. 

Results of our calculation are shown in Fig.l. The 
three dotted lines are the thermal photons from initial 
QGP at temperature 350,400 and 450 MeV respectively. 
In Fig.l, pQGD (hard) photons are shown as the 
dash-dotted line. Sum of thermal and pQCD photon 
yield are shown as the solid lines. pQCD photons alone 
can not explain the PHENIX data throughout the pr 



FIG. 1: PHENIX (preliminary) transverse momentum dis¬ 
tribution of direct photons in 0-20% centrality Au-|-Au colli¬ 
sions. The dotted lines are thermal photons from QGP with 
initial temperature 350,400 and 450 MeV (bottom to top), 
thermalised at ri=0.6 fm . The dash-dotted line is the pQCD 
photons. Solid lines are the sum of the pQCD and thermal 
photons. 


range. It explains the data at large pt but underesti¬ 
mate the yield at low pr- PHENIX (preliminary) di¬ 
rect photon data require thermal photons. Thermal pho¬ 
tons dominate the photon spectra at low px- As seen in 
Fig.l, if the initial QGP is formed at temperature ri=350 
MeV, the PHENIX data are underpredicted. The data 
are over predicted for initial temperature ri=450 MeV. 
PHENIX (preliminary) direct photon data are explained 
if the initial temperature is Ti=400 MeV. The analysis 
suggest that the PHENIX (preliminary) direct photon 
data in 0-20% centrality Au-|-Au collisions are explained 
if Au-|-Au collisions produce QGP with central tempera¬ 
ture 400 MeV at an initial time of 0.6 fm. 

Let us now consider the photon yield in the second 
scenario when Au-|-Au collisions produce hot hadronic 
gas. As it is well known, photon emission rate from 
QGP and from hot hadronic gas are very similar. De¬ 
generacy of resonance hadron gas is also of the same or¬ 
der as that of a QGP. Thus it is expected that a hot 
hadronic gas thermalised at ~ 0.6 at an initial tem¬ 
perature of Ti ^400 MeV will give a similar description 
to the PHENIX data, as it is obtained in the 1st sce¬ 
nario. However hadronic gas at such a high temperature 
is physically unacceptable. Density of the gas is very 
large ptadron ~ 50/m“^. Hadrons can not retain their 
identity at such a high density. What is the acceptable 
limit of initial temperature of the hot hadronic gas? Re¬ 
cent lattice QGD calculations indicate that critical tem¬ 
perature of confinement-deconfinement phase transition 
is Tc = 190 ± 10 MeV |2^. Possibly hadrons can re- 
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FIG. 2: PHENIX (preliminary) transverse momentum dis¬ 
tribution of direct photons in 0-20% centrality Au-|-Au colli¬ 
sions. The dotted lines are thermal photons from initial hot 
hadronic gas at initial temperature of 220 MeV, thermalised 
at ri=3,4 and 5 fm (from bottom to top). pQCD photons are 
shown as the dash-dotted The solid lines are the sum of the 
pQCD and thermal photons. Photons from a pre-equilibrium 
QGP are shown as the dash-dot-dot line. 


tain their identity still at higher temperature. We choose 
Ti=220 MeV as the physically acceptable value for the 
initial temperature of the hot hadronic gas. At this tem¬ 
perature density of the hadron gas is Phadrom 1.7fm~^ 
just below the limit phadron < 2fm~^ so that the hadrons 
retain their identity. Now what should be the thermalisa- 
tion time scale for a hot hadronic gas? Will it be as small 
as that of a QGP? If the hadronic gas at initial tempera¬ 
ture of 220 MeV thermalises at the same time scale (0.6 
fm) the PHENIX (preliminary) direct photon data are 
not explained. However, thermalisation time scale for a 
hadronic gas can be longer. Indeed, very small thermal¬ 
isation time scale (r^ ~ 0.6/m) obtained for a QGP is a 
puzzle in heavy ion physics and is not understood prop¬ 
erly. If the parton-parton collisions are responsible for 
thermalisation, the thermalisation time is significantly 
longer. Calculations performed within the ’bottom-up’ 
thermalisation scenario inclu ding 2 2 and 2 <-> 3 pro¬ 
cesses estimate > 2.6/m |^. Only mechanism 

for short thermalisation scale is plasma instability in soft 
modes of the gluon field [^. Due to rapid longitudinal 
expansion of the system, the momentum spectrum of par- 
tons quickly become anisotropic, the width of the longitu¬ 
dinal momentum distribution become narrower that the 
width of the transverse transverse momentum distribu¬ 
tion < Ap| ><<< >. A transverse chromo-electric 

field develop, which further enhances the fluctuations in 
the parton momentum distribution. Anisotropic parton 
momentum distribution can cause magnetic (transverse) 


instabilities, which can be very efficient in thermalising 
the system. Characteristic inverse time scale for insta¬ 
bility development is roughly of the order of gT for suf¬ 
ficiently anisotropic momentum distribution. Thermal¬ 
isation time scale of a hot hadronic gas is expected to 
be greater than that of a QGP as it is likely that in a 
hot hadronic gas, thermalisation will be driven by col- 
lisional process rather than by the plasma instabilities. 
Plasma instabilities, in hot hadronic gas, being electro¬ 
magnetic in nature, will not be efficient to thermalise the 
system rapidly (a << Og). In thermalisation of lin¬ 
ear sigma model fields were studied numerically. In the 
model, sigma model fields interact with a heat bath. Ir¬ 
respective of the initial field configuration, sigma model 
fields thermalises in the time scale ^ 5 fm. In resonance 
hadronic gas thermalisation process will be faster as num¬ 
ber of fields greatly exceed that of the sigma model. A 
reasonable estimate will be ~ 1-5 fm. 

In Fig.2 photon yield from the hot hadronic gas, with 
initial temperature Ti=220 MeV for three time scale of 
thermalisation, ri=3,4, and 5 fm are shown (the dotted 
lines). Even with large thermalisation scale, Ti=5 fm, 
PHENIX data are not explained if Au-|-Au collisions pro¬ 
duces hot hadronic gas at an initial temperature of 220 
MeV. The data remain underpredicted by a factor of two 
or more. The hadronic state can not produce the required 
number of photons. 

With large thermalisation time, =5 fm, the fluid 
matter spends considerable time in the pre-equilibrium 
stage. Thus while pre-equilibrium photons may not be 
important if Au-|-Au collisions produces QGP (rj=0.6 
fm), they may be important if the collisions lead to 
hot hadronic gas formation (Ti=5 fm). To unequivo¬ 
cally reject the hot hadronic gas scenario it is important 
to estimate the pre-equilibrium photons. Unfortunately 
emission of photons from a pre-equilibrium hadronic gas 
is not studied. However, photon emission from a pre¬ 
equilibrium QGP has been studied earlier . Pre¬ 

equilibrium photons are order of magnitude less than 
equilibrium photons. Following |29|| we have estimated 
the photon yield from a chemically non-equilibrated par- 
tonic system with initial conditions dictated by the HI- 
JING simulation for RHIG Au-|-Au collisions. The par- 
tonic system achieved kinetic equilibrium by the time 
Uso=0.31 fm at temperature of 570 MeV. At Tiso=0.31 
fm, gluon and quarks fugacities are 0.09 and 0.02 re¬ 
spectively. We have considered only longitudinal expan¬ 
sion. In Fig.2, pre-equilibrium photons (integrated over 
the time scale 0.31 to 5 fm) are shown as the dashed 
line. Photons from pre-equilibrium stage contribute in¬ 
significantly. The reason can be understood easily. The 
photon emission rate is weighted by the fugacities which 
remain at low values in the time scale integrated. With 
transverse expansion, pre-equilibrium emission at large 
Pt will increase, however, the increase would never be 
large enough to reckon them. If the pre-equilibrium hot 
hadronic gas contribute to the same order as the pre¬ 
equilibrium QGP, we can safely ignore them. Even if pre- 
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equilibrium photons from the hot hadronic gas exceed 
that from the pre-equilibrium QGP by a factor of 10, they 
still can be neglected. We conclude that PHENIX (pre¬ 
liminary) direct photon data are not explained if Au-|-Au 
collisions produces a physically acceptable hot hadronic 
gas, with initial temperature less than 220 MeV at an 
initial time less than 5 fm. 

To summarise, we have analysed the preliminary 
PHENIX data on the transverse momentum distribution 
of direct photons in 0-20% centrality Au-|-Au collisions 
at y^SNN=200 GeV. Two scenarios are considered, (i) 
Au-I-Au collisions produces a QGP and (ii) Au-|-Au colli¬ 


sions produces hot hadronic gas. PHENIX (preliminary) 
direct photon data are explained in the first scenario if 
QGP is produced at initial time Ti=0.6 fm at an ini¬ 
tial temperature ri=400 MeV. PHENIX data are not 
explained in the alternate scenario if the hot hadronic 
gas is produced with initial temperature less or equal to 
220 MeV at an initial time 5 fm or less. The data remain 
underpredicted. As the hot hadronic gas at higher tem¬ 
perature or with longer thermalisation time is physically 
unacceptable, we conclude that PHENIX (preliminary) 
direct photon data are explained only if QGP is produced 
in 0-20% centrality Au-|-Au collisions. 
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